function rest=addconv(alpha,m)

J=indexJ(k,n);
numofindex=size(J,1);
K=size(J,2);
totalequation = numofindex*(2*m+1);
beta=zeros(1,numofindex);

rest=0;
for pp=1:numofindex
        p=indexJ(pp,:);
        betas=getbeta(alpha);
        for bb=1:size(betas,1)
            beta=betas(jj,:);
            u1n=findindex(J,beta+p);u2n=findindex(J,alpha-beta+p);
            if(u1n&&u2n)                
                rest=rest+coefficientC(alpha,beta,p) ...
                    *convlution(u((u1n-1)*(2*m+1):u1n*(2*m+1)),...
                                u((u2n-1)*(2*m+1):u2n*(2*m+1)),m);
            end
        end
end
rest=-1i*pi*(-m:m)'.*rest;